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Abstract 

The square lattice spin-1/2 antiferromagnet containing a dilute concentration, 5, of ran- 
domly placed ferromagnetic nearest-neighbor bonds is studied at low-temperature via non- 
linear cT-model techniques and by exact diagonalization. We generally find that long-range 
Neel order is destroyed above a critical strength in the defective ferromagnetic exchange 
coupling-constant given by \Kc\/J ~ We also observe large statistical fiuctuations 

both in the spin-stiffness and in the antiferromagnetic structure-factor near this critical 
point, suggesting the onset of a spin-glass phase. 
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It is now well established that the spin-1/2 Heisenberg model on the square-lattice 
with nearest-neighbor antiferromagnetic coupling displays long-range Neel order at zero- 
temperature, along with a correlation length that diverges exponentially at low-temperature.^ 
This model, in fact, successfully accounts for the long-range two-dimensional (2D) anti- 
ferromagnetic correlations that are observed in "parent" compounds to high-temperature 
superconductors such as La2Cu04 at temperatures just above the transition temperature 
to three-dimensional Neel ordering. -"^ However, the latter antiferromagnetic insulating state 
disappears rapidly upon hole doping, where it is followed eventually by the appearance of 
the superconducting phase. ^ 

Aharony et al. have proposed that such doping results in a frustrating eflFective fer- 
romagnetic coupling between the nearest-neighbor copper spins situated near the hole, 
which is presumed to lie statically at the intermediate oxygen super-exchange site.^ (A 
three-band Hubbard model is implicit in this approach.) The observation of a disordered 
spin-state via neutron scattering experiments in the La2-a;Sra;Cu04 series for temperatures 
below 100 K, and in the doping range 0.015 < a; < 0.05, supports this scenario.^ In this 
Letter, we elucidate the theoretical consequences of the former proposal by studying the 
low-temperature properties of a spin-1/2 nearest-neighbor Heisenberg antiferromagnet on 
the square-lattice with a quenched random distribution of dilute ferromagnetic nearest- 
neighbor bonds. In particular, we first consider the corresponding classical problem at 
non-zero temperature via the non-linear u-model,'^ which is known to correctly describe 
the pure 2D antiferromagnet.^ We then compute the spin-stiffness, as well as various spin- 
correlations, of the full quantum problem at zero-temperature by exact diagonalization 
on a 4 X 4 lattice using the Lanczos technique.^ For a fixed concentration 5 of frustrating 
ferromagnetic bonds, both studies find that long-range Neel order disappears on average at 
a critical strength in the nearest-neighbor ferromagnetic exchange coupling-constant given 
by \Kc\/J ~ 5~^/'^. In the case of antiferromagnetic defective bonds, however, the exact 
diagonalization study finds no evidence for the destruction of long-range Neel order, as one 
intuitively expects. In addition, the statistical fiuctuations shown by both the spin-stiffness 
and the antiferromagnetic structure-factor obtained in the exact diagonalization study be- 
come extremely large near the former critical point, suggesting the onset of a spin-glass 
phase. ^ Last, our numerical calculations also reveal that uniform chiral spin-fluctuations^ 



are enhanced by the introduction of either antiferromagnetic or ferromagnetic defective 
bonds. A discussion of these results within the context of the low-temperature magnetic 
phase diagram of La2-xCu04 is given at the end of the paper. 

Consider a nearest- neighbor spin- 1/2 Heisenberg model Hamiltonian, 

H = ^JijSi - Sj, (1) 

(y) 

on the square-lattice, with a dilute concentration, S, of randomly distributed bonds, Jy = 
K < 0, lying within a uniform antiferromagnetic background with Jij = J > OJ Clearly 
then, the average value and the standard deviation of the exchange coupling-constant per 
bond are given by 

Jo = KS + J{l-5) and Ji = \J-K\[{l-S)S]^/^, (2a, 6) 

respectively. Below, we analyze the classical non-linear cr-model for the present antiferro- 
magnet with random frustration and find that the destruction of long-range Neel order at 
low temperature on average is given by the Lindemann-like criterion Jq ~ J^. 

Classical Non-linear a -model. Suppose that long-range Neel order exists in the ground 
state of Hamiltonian (1). In analogy with the treatment of the pure 2D quantum antiferromagnet,^| 
we consider then the classical 2D non-linear cr-model^ with a randomly fluctuating local 
spin-stiffness field, which presumably mimics the random frustration present in Hamilto- 
nian (1). The partition function for this model may be expressed as 



Z = j Vn{x)VX{x')e^p{-pE), (3) 
with an energy functional given by 

E=Y j d''^{[Jo + Mm^n{x)\^+i\{x)[\n{x)\^-l]} (4) 

such that Ji (x) is a randomly fluctuating local spin-stiffness fleld with zero average- value. 
Here, fi{x) is the Neel order parameter constrained to unit modulus by integration over 
the Lagrange-multiplier fleld A(x), while s = 1/2 is the spin per site. (Quantum effects 
renormalize down in the pure case.^'^'^) Fourier transformation of all the fields above 
yields a fluctuation contribution to the energy (4) equal to 

^1 = - y E E E MQ)[m ■ n{k')]{k ■ k')S{k + k' + q). (5) 



For the sake of calculational convenience, let us now generalize the three-dimensional Neel 
order-parameter to an A'"-dimensional vector field. In the large-A?" limit we may ignore 
fluctuations in the Lagrange-multiplier field and take X{x) = Aq.*'^ If the effect of fiuc- 
tuation energy (5) is accounted for within the coherent potential approximation (CPA),^ 
which ignores backscattering, then the Greens function for the Neel order-parameter field 
is given by {na{k)np{-k)) = G{k)Sai3, with 

G{k) = {PsY^[Jok^ + iXo - (6) 

where the self-energy correction is self-consistently determined by 

m = Ips'a'jfJ -^{k.k'fGik'). (7) 

Above, it has been implicitly assumed that {Ji{x)Ji{x')) = 2oF'Jl5{x — x'), where the 
constant a corresponds to the lattice-spacing of the original Heisenberg model (1). In 
general, the constraint that the Neel order parameter be of unit modulus must also be 
satisfied;^ i.e., we have that 1 = (n^) = N J{2TT)~'^(fkG(k). Substitution of the ansatz 
S(fc) = 7 JqA;^ into both this constraint and Eq. (7) yields the relationships 

(1 - 7)X = ln{kl^'^ + 1) and (8a) 

(1 - = Y[i - {kncrMkie + 1)] m 

respectively, where the antiferromagnetic correlation length, ^, is defined by iXo = (1 — 
7)Jo^~^, where X = Atts^Jq/NT, and where Y = {kj^a'^ /47T){jf /4J^). Here, ko rep- 
resents the momentum integration cut-off that we henceforth take to satisfy k'j-,a? /An — 
{2n)~'^ d^k = 1, by correspondence with the first Brillouin-zone of the square-lattice. 
To begin the analysis of the above mean-field Eqs. we consider first the low-temperature| 
regime, such that X ^ 1. Then Eq. (8a) implies that the correlation length diverges ex- 
ponentially, like in the pure system, ^'^ if 7 < 1. Hence, by Eq. (8b) we have that 

(1 -^)^^Y (9) 

in this case. Yet in the dilute limit, 5 — > 0, we also have that y <S 1 by Eq. (2b). Hence, 
(9) implies that 7 = y. This result in conjunction with (8a) yields, therefore, that the 
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correlation length is approximately given by k^$,^ = e^^~^^'^ at low temperature. Within 
CPA, therefore, the effect of dilute random fluctuations in the local spin-stiffness is to 
reduce the effective spin-stiffness connected with long-range Neel-order-'^''^ by a factor of 
1 — Y — 1 — JI/AJq. On the other hand, (9) also yields that the parameter Y has a 
maximum value of 1/4 at 7c = 1/2. The CPA approximation thus also results in a critical 
ferromagnetic coupling-constant, Kc < 0, below which solutions to mean-field Eqs. (8a) 
and (8b) possessing long-range Neel-order no longer exist, determined by the condition 
Jq = . Substitution of expressions (2a,b) into this condition then gives 

li^clcPA = J/S'/^ (10) 

in the dilute limit. The exact solution for the critical value of the ferromagnetic exchange 
coupling-constant within CPA is displayed in Fig. 1 as a function of concentration. It is 
important to remark that the spin-stiffness precisely at this critical point does not vanish, 
but is reduced by a factor of 1 — 7c = 1/2 with respect to stiffness of the pure system. 
Second, a straight-forward analysis of mean-field Eqs. (8a) and (8b) demonstrate that the 
paramagnetic solution ^ — > is possible either at high-temperatures or when the spin- 
stiffness vanishes. In particular, one finds that 7 and — > (1 — \XY)X <^ 1, 
which implies either that X <^ 1 or that 7 \XY — > 1. Last, it can also be shown that 
no zero-temperature solution to mean-field Eqs. (8a,b) exists that has a non-zero yet finite 
correlation length. 

Exact Diagonalization. In ref. 5, the authors have implemented a direct method for 
calculating the spin-stiffness of generic spin-1/2 Heisenberg models (1) by exact diagonal- 
ization of small clusters with spin-dependent twisted boundary conditions in the = 
subspace.-^° The method determines the spin-stiffness of the ferromagnet to an accuracy of 
one percent, while it is accurate to within ten percent in the case of the antiferromagnetic 
chain. ^ We have applied this technique to the present Hamiltonian (1) describing dilute 
randomly placed defective bonds within a background of homogeneous antiferromagnetic 
bonds on a 4 X 4 square-lattice. Fig. 2 displays both the average and the standard devia- 
tion of the so-determined spin-stiffness for 100 randomly chosen configurations of A*";, = 2, 4 
defective bonds as a function of the coupling-constant, K. The first observation to be made 
is that the average spin-stiffness decreases rapidly to zero with increasing ferromagnetic 
coupling, ■^■'^ \K\, while it remains positive if the defective bonds are antiferromagnetic. The 
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critical value of K at which this average value vanishes is plotted in Fig. 1 as circles for 
the case of Ni, = 1, 2, 3, 4 and 5 defective bonds. The agreement with the previous classical 
non-linear cr- model results is quite good, indicating that quantum effects do not play an 
important role in the destruction of long-range Neel order by such random frustration. 

The next important feature to observe is that the standard deviation of the spin- 
stiffness obtained via the Lanczos study, in contrast to its average value, increases dramat- 
ically with increasing strength of the ferromagnetic bonds (see the inset to Fig. 2). The 
presence of such large statistical fluctuations in the spin-stiffness supports the claim that a 
spin- glass phase appears once long-range Neel order is destroyed, as was first suggested by 
Aharony et al..^ Such a state is also consistent with our findings that the ground-state ob- 
tained in the present exact-diagonalization study is a spin-singlet in the transition regime, 
while the first excited state in this regime has spin-one. Last, we mention that recent 
theoretical studies of quantum Heisenberg models (1) with random infinite-range defective 
bonds also find evidence for spin-glass behavior. 

To complete our numerical investigation, we have measured the spin-spin correlation 
function (Si-Sj) for the ground-state of Hamiltonian (1) obtained by the Lanczos technique. 
Fig. 3 displays both the average and the standard deviation of its Fourier transform, S{q), 
at the antiferromagnetic point q = (tt, tt) for 100 randomly chosen configurations of sys- 
tems with Nb = 2,4 defective bonds as a function of the coupling-constant, K. Contrary to 
the spin-stiffness (Fig. 2), we see that both ferromagnetic and antiferromagnetic defective 
bonds suppress long-range antiferromagnetic spin correlations. The fact that <S'(7r, tt) re- 
mains positive near the critical point where the spin-stiffness vanishes on average could be a 
finite-size effect. Fig. 3 also shows that the standard deviation of S{7r, tt), which is a type 
of Edwards- Anderson spin-glass order-parameter,-^^ increases dramatically upon the intro- 
duction of frustration, eventually saturating to a value on the order of the structure-factor 
itself. If we presume that the spin-glass onset point coincides with the half-maximum of 
this standard deviation, then we obtain critical values for the defective exchange coupling- 
constant of \Kc\ ~ 2,4 in the cases of four and two frustrating bonds, respectively. This 
is consistent with the point at which long-range Neel order is destroyed, as indicated in 
Fig. 1, and with the proposal that spin-glass order follows.^ We have also computed the 
average of the squared uniform chiral moment, which we define on an arbitrary plaquette 
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of points z = 1,2,3,4 by = [Si ■ {S2 - S4) x + S2 ■ {S3 - Si) x 84]'^. A straight 
forward yet tedious calculation reveals the remarkable identity that = i(l — P13P24:), 
where Pij denotes the spin-exchange operator between sites i and j. (Notice that the Neel 
configuration is a null eigcnstate of this operator!) Our results for this chiral moment are 
plotted in the inset to Fig. 3. In contrast to the behavior shown by the antiferromagnetic 
correlations S'(7r,7r), we observe that generally increases with increasing strength in 
both antiferromagnetic and ferromagnetic defective bonds. It is worth mentioning here 
that staggered chiral spin-fluctuations are exponentially suppressed at low temperature in 
the pure quantum-renormalized 2D Neel phase.® 

In conclusion, both the non-linear cr-model and direct calculations of the spin-stiffness 
for 2D spin-1/2 Heisenberg models based on exact-diagonalization of a 4 x 4 lattice with 
twisted boundary conditions^' show evidence for the rapid destruction of long-range Neel 
order at low-temperature induced by the presence of dilute ferromagnetic bonds. In ad- 
dition, the appearance of large statistical fluctuations in the latter spin-stiffness near the 
point where its average value vanishes support the claim that the disordered phase is a 
spin-glass.^ Hence, the low-temperature phase diagram of our Hamiltonian (1) appears to 
be quite different from that of the pure quantum non- linear cr-model in 2 -|- 1 dimensions,-^ 
which has recently been discussed in connection with the high-temperature magnetic be- 
havior displayed by the doped "parent" compounds to high-Tc superconductors.^^ Last, 
the low-temperature Neel phase in La2-a;Cu04 disappears at a hole concentration of 
X = 0.015.^ If we suppose that each of these holes corresponds to a ferromagnetic bond 
in the present model (1),^ then a critical concentration 5 = x/2 = 0.0075 of such bonds 
implies a ferromagnetic exchange coupling-constant of |K| ~ 10 J by the previous CPA 
result (10). The latter value is consistent with simple microscopic estimates based on di- 
rect copper-oxygen exchange in the oxide superconductors.^ Furthermore, inversion of Eq. 
(10) shows that the critical concentration for the destruction of long-range Neel order is 
very sensitive to \K\/J, suggesting that the low-temperature magnetic phase diagram of 
the oxide superconductors could be sensitive to pressure. Last, we emphasize that these 
results indicate that long-range Neel order exists up to a non-zero critical concentration of 
frustrating bonds, contrary to statements in the literature suggesting that any amount of 
random frustration is sufficient to destroy Neel order. ^ 
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Figure Captions 



Fig. 1. The solid line given by the solution to the quadratic equation Jq = [see Eqs. (2a,b)] 
delineates the point beyond which long-range Neel order is no longer possible in the 
present cr-model calculation. The circles correspond to the point at which the spin- 
stiffness vanishes on average in the 4x4 exact diagonalization study for A^^ = 1, 2, 5 
ferromagnetic bonds, with the concentration determined by S = Ni,/32. 

Fig. 2. Shown is the average spin-stiffness for A^^ = 2, 4 defective bonds as a function of the 
defect exchange coupling-constant K. The standard deviation of the spin-stiffness is 
displayed in the inset. 

Fig. 3. The average long-range antiferromagnetic spin correlations are plotted as a function of 
the defect coupling-constant K for Ni, = 2,4 defective bonds, while the corresponding 
plot for the squared uniform chiral moment is shown in the inset. Error bars represent 
the standard deviation in each quantity. 
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